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Abstract 



In quantum adiabatic evolution algorithms, the quantum computer follows the 
ground state of a slowly varying Hamiltonian. The ground state of the initial 
Hamiltonian is easy to construct; the ground state of the final Hamiltonian 
encodes the solution of the computational problem. These algorithms have 
generally been studied in the case where the "straight line" path from initial 
to final Hamiltonian is taken. But there is no reason not to try paths involving 
terms that are not linear combinations of the initial and final Hamiltonians. 
We give several proposals for randomly generating new paths. Using one of 
these proposals, we convert an algorithmic failure into a success. 

1 Introduction 

& . 

Quantum adiabatic evolution algorithms |J are designed to keep the quantum com- 
puter in the ground state of a slowly varying Hamiltonian H(t). The initial Hamil- 
tonian, Hp = H(0), is chosen so that its ground state is easy to construct. The 
ground state of the final Hamiltonian Hp = H(T), where T is the running time of 
the algorithm, encodes the solution to the computational problem at hand. A simple 
way to construct the interpolating Hamiltonian is to define 

H(t) = H(t/T) < t < T (1) 

where 

H(s) = (1 - s)H B + sH P < s < 1 . (2) 

But there is no reason to restrict attention to this path from Hp, to Hp. The pos- 
sibility of varying the path was discussed in [Q, 0, and M. The adiabatic algorithm 
will work taking any path H(s) with H(0) = Hp and H(l) = Hp as long as T is 
much larger than 1 / gap 2 where the gap is the minimum energy difference, as s varies, 
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between the ground and first excited states of H(s). In this paper we consider paths 
of the form 



We view He as an extra piece of the Hamiltonian that is turned off at the beginning 
and end of the evolution. For a given instance of a problem, that is, for a fixed Hp, 
we imagine running the algorithm repeatedly with different choices of H E , possibly 
chosen randomly. 

In this paper we propose some choices for H E . We then apply one of these choices 
to a problem discussed in || and show that the addition of He can change the 
performance of the algorithm from unsuccessful to successful. 

2 Examples of different paths 

We are most optimistic about the performance of quantum adiabatic algorithms when 
they are applied to the classical problem of finding the minimum of a "local" cost 
function. By this we mean a function h(zi, ... ,z n ) where each Zi = 0, 1 and 



and each he is a nonnegative integer- valued function that depends only on a few Zj. 

The Hamiltonian that governs the quantum algorithm is constructed as follows. 
For each term he we define an associated quantum operator Hpc- Suppose, for 
example, that he depends on three bits z'c, jc, and fee- Then we have 




(3) 




(4) 



c 





Accordingly we define the problem Hamiltonian 




(6) 



c 



For the beginning Hamiltonian we define 





where (again) i c , jc, and k c are the three bits appearing in h c . Now 




(8) 



c 



and using @ and (0) we see that H(s) from (0) can be written 




(9) 



c 
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This is the form of H(s) that was originally proposed in equations (2.26) and (2.28) 
in 0. 

We now turn to H E . To preserve the J^c on ^ ne right-hand side of (|9|), we pick a 
term H E q for each term Hq. This H E q will only depend on the bits involved in he- 
We then have 

#E = ^#E,C (10) 

c 

and H(s) from (|3|) is now 

H(s) = J^ii 1 ~ s ) H b,c + sH PtC + s(l - s)# E ,c] . (11) 
c 

Note that for each clause the coefficients (1 — s), s, and s(l — s) could all be replaced 
by any three smooth functions with the same corresponding values at s = and s — 1, 
but we make the clause-independent choices (1 — s), s, and s(l — s) for simplicity. 
Finally we must choose H E q. Here are three proposals. 

PI. For each term he involving 6c bits let Aq be a randomly chosen Hermitian 
matrix of size 2 bc x 2 bc . We imagine that the distribution of Aq produces 
entries with magnitude of order unity. We also propose setting the diagonal 
entries in each Aq to 0. Now i?E,c is defined by letting the matrix Aq act 
on the be bits involved in he- The diagonal avoids confusing the effects of 
varying the path in the quantum algorithm with modifications of the classical 
cost function itself. See ||]. 

P2. For this proposal we specialize to the case where each term he has the same 
form independent of C. For the case of three bits this means that there exists 
an /i 3 such that 

hc( z ici z jci Z kc) = hz\ Z ici Z 3Ci Z kc) (1^) 

for each C. (An example is the NP-complete problem Exact Cover where 
h 3 (z,z',z") = 1 if z + z' + z" ^ 1 and if z + z' + z" = 1.) For simplicity 
we restrict to the case where each he depends on three bits. We propose gen- 
erating a single (order unity) random 8x8 Hermitian matrix A and then using 
this for each i?E,c- Each H E q has the 8x8 matrix A acting on bits ic, jc, an d 
kc- Again, we propose choosing the diagonal entries of A to be 0. 

P3. Here we make a proposal for 3-SAT. An n-bit instance of 3-SAT is a collection 
of clauses. Each clause specifies which three bits are in the clause. Each clause 
is True for seven of the eight assignments of the three bits and False on one 
assignment, so there are eight different types of clauses. Accordingly, he = 1 if 
(zj c , Zj c , Zk c ) is the False assignment and he = on the other seven assignments. 
Thus h(zi, . . . , z n ) is the number of clauses violated by the string z%, . . . ,z n and 
a satisfying assignment exists iff the minimum value of h is 0. 
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We can view the eight types of 3-SAT clauses as variants of a single basic clause, 
say the one for which (0, 0, 0) is the False assignment. The other seven clauses 
are related to the basic clause by negation of some (or all) of the bits. For this 
proposal for we start with a single 8x8 matrix A and associate it with the 
basic clause. For the other seven clauses we have seven other 8x8 matrices. 
Each is obtained from A by negating the same bits that are negated in the 
corresponding clause. (Bit negation of A is achieved by conjugating A with the 
o x matrix acting on the bit.) A might be chosen randomly but again we take 
the diagonal entries to be 0. This proposal has the invariance property that if 
an instance of 3-SAT is changed by replacing some bits by their negations in all 
clauses, the performance of the quantum algorithm is identical. 



3 A different path can turn failure into success 



Here we reexamine an example given in || where the quantum adiabatic evolution 
algorithm, evolving according to (|2|), without an ifg, was shown to require a time 
exponential in the number of bits. We will now show that for chosen randomly 
according to Proposal P|| above, the quantum adiabatic algorithm succeeds in poly- 
nomial time. 

The example has a cost function h, of the form (|4]), where each he depends on 
three bits and is of the same form for each C, that is, equation ([12]) holds. In this 
example 



h 3 (z,z',z") 



( z + z' + z" = 

3 z + z' + z" = 1 

1 z + z' + z" = 2 

1 z + z' + z" = 3 



(13) 



The cost function includes a term for each set of three bits giving 



h(z 1 , ... ,z n ) = ^ h 3 { 

i<j<k 



(14) 



Note that the global minimum of h occurs at z% — 0, Zi — 0, . . . , z n — where h is 0. 

The performance of the quantum adiabatic algorithm on this highly symmetrized 
problem can be analyzed asymptotically in the number of bits n. Using the symmetry 
and equations ®, (H), (|i~3D , and (|14|) we can write 



-{--s z 

2 V2 



1 (n 
+ 2U 



1 



+ -(--s t 

6 V2 



»_ S ,_l)(»_ S ,_2) (15) 
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where S z is the ^-component of the total spin: 

n 
i=l 

For H-q we use (|7]) and (^) to obtain 

2 yv2 



^H^f?-^) (IT) 



where 



(0 

i=l 



and for later use 



5 w = |£aW. (19) 



We now turn to i?E generated by Proposal P|2|, from a single 8x8 matrix A with 
0's on the diagonal. By ( JlQ| ) and the fact that all sets of three bits are included in the 
sum we find that is a cubic function of S x , S y , and S z whose coefficients are linear 
combinations of the entries of A. The choice that A is on the diagonal implies that 
there are no S z , S^, and S% terms in whereas these are the only terms in Hp. 

In this symmetric case the total Hamiltonian H(s) given by (O) is a function of 
the total spin operators S x , S y , and S z . For n large, we make the ansatz that for 
each s the ground state of H(s) is a state W, </?) that is an eigenstate of the total spin 
in the (9, <p) direction 

(sin 9 costp S x + sin 9 sin ip S y + cos 9 S z ) \ 9,Lp^ = —\0,ip) ■ (20) 

With this ansatz finding the ground state of H (s) is reduced to finding the minimum 
over (9, ip) of 

(9,<p\H(s)\9,<p) . (21) 
To evaluate (BTJ) in the large-n limit we use relations such as 
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(0,(f\S x 


9,<p) = 


n 

— sin 

2 


{9,if\Sy 


9,<f) = 


n 

-sin 


(9 co\S 


9 in) = 


n 

— cos 


(0,<p\S£ 


o,v) = 




(9,<p\S X Sy 


e,v>) = 




(9,cp\S x S z 


9,<p) = 


(I) 2 



/Tl\ 

[-) sin 2 9 cos 2 f + 0(n) 

sin 2 # cos 93 sin 99 + 0(n) 
sin 9 cos cos 90 + O (n) 



(22) 



We now define the effective potential 



2\3 



V(9,<p,8)= lim(-) {0,<p\H{s)\e,<p) 



which we break into two parts: 



V (9, if, s) = lim (-) 3 (9, <p\{(l- s)H B + sH P }\9, if) 



and 



V E (9,if,s) = lim f-) 3 (9,if\s(l-s)H E \9,if) 



(23) 



(24) 



(25) 



so that 

V = V + V E . (26) 
We now review || what happens without H E . We have 

V {9,if,s) = 2(1 - s)(l - sin# cosf) + ±s(13 + 3 cosfl - 9 cos 2 9 - 7 cos 3 0) . (27) 

The minimum of Vo for any s is at if = 0. In Figure |1] we plot Vo(#, 0, s) for six values 
of s. In the first panel s = and the effective potential is minimized at 9 = 7r/2, which 
corresponds to the ground state of i^B- As s increases the minimum moves to the right 
(see s = .2). Continuously following this local minimum leads to the local minimum 
at 9 = 71 when s — 1. A quantum computer running the adiabatic algorithm for a 
time only polynomial in n will evolve close to the state \9 — 7r, if — 0) at the end of 
the run. This state \9 = 7r, y> = 0) = |zi = 1> ^2 = 1> ■ ■ ■ , z n = l) does not correspond 
to the ground state of Hp, which is at \9 = 0, ip = 0) = \Zi = 0, z% = 0, . . . , z n = 0). 

The failure of the quantum adiabatic algorithm on this highly symmetrical prob- 
lem with n 3 clauses may not be a good test of the effectiveness of the algorithm on 
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s=0 



s=.2 



s=.4 






Figure 1: The diamond continuously tracks a local minimum of the effective potential, 
corresponding to the behavior of the quantum adiabatic algorithm run for polynomial time. 
The final point reached is not the global minimum of the effective potential at s = 1 and 
the algorithm fails. Here He is absent. 



computationally interesting problems. But in any event this failure can be turned to 
success with the addition of H^. 

Although Proposal PfJ uses a random 8x8 matrix A, we begin by studying the 
effects of i?E arising from a carefully chosen A: 



A 



which gives rise to 
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(28) 



H E = -2n(S x S z + S z S x ) + 0{n 2 ) 



(29) 
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from which we get 



Ve(8, <f, s) = — 8s(l — s) cos 9 sin 9 cos <p 



(30) 



The full effective potential V(9, (p, s) given by ( |2"o| ) with Q27D and ( 50 ) can only be 
extremized at (p = and 7T. In Figure |2] we show V(f9, 0, s) as a solid line and 
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Figure 2: The diamond continuously tracks a local minimum of the effective potential for 
the Hamiltonian including a particular He- (The solid line corresponds to p> = 0; the dotted 
line to ip = it.) The final point reached is the global minimum of the effective potential at 
s = 1, and the algorithm succeeds in polynomial time. 

V(0,ir,s) as a dashed line for six values of s. At s = the ground state of H(0) is 
\d = 7r/2, (p = 0). In the first panel we see that # = 7r/2,(p = 0is the global minimum 
of V(9, (p, 0). As s increases, the global minimum moves continuously to 9 — 0, (p — 0. 
The state = 0, <p = 0) = \z\ = 0, Z2 = 0, . . . , z n = 0) is the ground state of Hp and 
corresponds to the minimum of the cost function given by (O) and (0). Running 
the quantum adiabatic evolution algorithm for a time polynomial in n will bring the 
quantum computer close to the state \9 = 0,cp = 0), solving the problem of finding 
the minimum of h. 

The example just given is not an isolated success. Clearly the same behavior will 
be seen for A's close to fl2"B|). Moreover, we performed the following experiment. We 



Quantum Adiabatic Evolution Algorithms with Different Paths 



9 



toss a random real symmetric 8x8 matrix A with nondiagonal elements independently 
and uniformly distributed between —3 and 3 and with the diagonal set to 0. We 
use A to construct H E and the corresponding Ve- We then continuously track a local 
minimum of V(9, if, s) as s varies from to 1, starting at 9 = n/2,ip = when s = 0. 
In 351 out of 1000 tries, this local minimum moves to 9 = 0, tp = 0, which is the 
global minimum of the effective potential at s = 1. Thus, in these tries, the adiabatic 
algorithm run for polynomial time will end up close to the ground state of Hp. This 
shows that a large "volume" of the space of possible A's leads to an He that ensures 
finding the global minimum of h in polynomial time. 

(It is conceivable that in some of these tries the continuously tracked local mini- 
mum, which coincides with the global minimum at s = and s — 1, does not coincide 
with the global minimum for all intermediate s. If this happens, the algorithm run 
for polynomial time will still succeed in finding the global minimum of h, because the 
probability of leaving the state corresponding to the local minimum is exponentially 
small.) 

Our analysis, via the effective potential, of the quantum algorithm with an He 
chosen according to Proposal P2 gives the asymptotic (large n) behavior of the al- 
gorithm. This analysis is only possible because of the symmetry: the same terms 
are included in the Hamiltonian for each set of three bits. The function h% in ( p~3| ) 
could also be viewed as arising from a sum of terms corresponding to thirteen 3-SAT 
clauses. We could then use Proposal P3, which is designed for 3-SAT. Here one could 
apply the effective potential analysis and determine the asymptotic behavior of the 
algorithm with random choices of A. Presumably the behavior is similar to what we 
found using proposal P2, but we have not done the analysis. 

Proposal PI could be used instead in this problem, but choosing a different Ac 
for each C breaks the symmetry and we don't know how to determine the large n 
performance of the algorithm in this case. In [Q a study of random 9-bit 3-SAT in- 
stances was performed using proposal PI. The gaps were compared with and without 
He- Overall the gaps with He were smaller than without He- However, the gaps that 
were smallest without He often increased after the addition of He and were no longer 
small compared to the other gaps. This argues for running the quantum adiabatic 
algorithm on each instance with a variety of randomly chosen paths. 



4 Conclusion 

Previous work on quantum adiabatic evolution algorithms has assumed a "straight 
line" path in Hamiltonian space from the beginning Hamiltonian Hp to the problem 
Hamiltonian H-p as in @. This path is simple but arbitrary. The idea of trying to 
remain close to the ground state of a time-varying Hamiltonian is equally applicable 
to other paths from Hp, to Hp, such as (j3|). The quantum adiabatic algorithm is 
guaranteed to give a 100% probability of success only as the run time goes to infinity, 
so a finite run time requires repetition. We suggest that repetitions be carried out with 
different paths. In cases where the gap in H(s) given by (£3) happens to be small, the 
required run time is correspondingly long. Adding He as in (Q) may avoid the region 
in Hamiltonian space where the gap is small and be computationally advantageous. 



10 



E. Farhi, J. Goldstone, S. Gutmann 



Acknowledgments 

We thank Michael Sipser and Wim van Dam for useful discussions. This work was 
supported in part by the Department of Energy under cooperative agreement DE- 
FC02-94ER40818 and by the National Security Agency (NSA) and Advanced Re- 
search and Development Activity (ARDA) under Army Research Office (ARO) con- 
tract DAAD19-01- 1-0656. 



References 

[1] E. Farhi, J. Goldstone, S. Gutmann, M. Sipser, "Quantum Computation by 
Adiabatic Evolution" , |quant-ph/0001 106j . 

[2] E. Farhi, J. Goldstone, S. Gutmann, J. Lapan, A. Lundgren, Daniel Preda, "A 
Quantum Adiabatic Evolution Algorithm Applied to Random Instances of an 
NP-Complete Problem" , |quant-ph/0104129| . 

[3] A.M. Childs, E. Farhi, J. Preskill, "Robustness of Adiabatic Quantum Compu- 
tation" , |quant-ph/010804Sl Phys. Rev. A65, 012322 (2002). 

[4] D. Preda, "Quantum Adiabatic Evolution with a Different Path". MIT Senior 
Thesis, June 2002. 



[5] E. Farhi, J. Goldstone, S. Gutmann, "Quantum Adiabatic Evolution Algorithms 
versus Simulated Annealing", quant-ph/ 020103"! . 



